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Understanding the origin of the orbital resonances of the Galilean 
satellites of Jupiter will constrain the longevity of the extensive 
volcanism on Io, may explain a liquid ocean on Europa, and may 
guide studies of the dissipative properties of stars and Jupiter-like 
planets. The differential migration of the newly formed Galilean 
satellites due to interactions with a circumjovian disk can lead to 
the primordial formation of the Laplace relation n x — 3n 2 + 2n 3 = 0, 
where the are the mean orbital angular velocities of Io, Europa, 
and Ganymede, respectively. This contrasts with the formation of 
the resonances by differential expansion of the orbits from tidal 
torques from Jupiter. 

The orbital resonances among the Galilean satellites of Jupiter have led to sustained 
dissipation of tidal energy to produce astounding volcanos on Io and, probably, to maintain 
a liquid ocean on Europa. Understanding the origin of these resonances will constrain the 
history of the satellites and their formation scenarios. The two models proposed are an 
assembly of the resonances through differential tidal expansion of the orbits from tides 
raised on Jupiter ([?],[§) and a primordial origin of unspecified assembly and subsequent 
evolution away from exact resonance Here we discuss a means of assembling the 

resonances during the formation process in support of a primordial origin. 

The currently observed orbital resonances at the 2:1 mean motion commensurabilities 
involving Io-Europa and Europa-Ganymede are such that the resonance variables Q\ = 
Ai — 2A 2 + w\ and #3 = A 2 — 2A 3 + ro 2 librate about 0° and 9 2 = \\ — 2A 2 + to 2 librates 
about 180°, all with small amplitude. Here the Aj are mean orbital longitudes and the 
Wi are the longitudes of periapse numbered consecutively for Io, Europa, and Ganymede, 
respectively. The simultaneous libration of 62 and #3 means the Laplace angle #5 = 
O2—Q3 = Ai— 3A 2 +2A3 ~ 180° with small libration amplitude, leading to n\— 3n 2 +2n3 = 0, 
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where rii = A» are the mean angular velocities with the dot indicating time differentiation. 
This last libration condition is called the Laplace relation because Laplace first understood 
its stability. The geometry is such that conjunctions of Io and Europa occur when Io is 
near the periapse of its orbit and Europa is near its apoapse, and conjunctions of Europa 
and Ganymede occur when Europa is near its periapse. Ganymede can be anywhere in 
its orbit at the conjunctions with Europa, because the resonance variable # 4 = A 2 — 2A 3 + 
©3 circulates through all angles. The existence of this group of resonances with small 
amplitudes of libration implies an assembly of originally random orbits by dissipative 
processes. 

A currently popular means of resonance assembly is by differential orbital expansion 
due to gravitational tides raised on Jupiter, principally by Io, and damping of the libration 
amplitudes by tidal dissipation within the satellites, mostly in Io Io's orbit expands 

more rapidly than Europa's, and capture of Io and Europa into the 2:1 orbital resonances 
is certain for a reasonably wide range of initial orbital eccentricities. Continued orbital 
expansion tends to drive the pair deeper into the resonance (i.e., driving n\ — 2n 2 to 
smaller values) and thereby to increase the orbital eccentricities. But dissipation of tidal 
energy within the satellites tends to decrease the eccentricities, so an equilibrium can be 
reached where the eccentricities and the ratio of the mean motions are almost constant 
as the pair of satellites moves out together locked in the resonances ([J,[|). Europa even- 
tually encounters the 2:1 mean motion commensurability with Ganymede, where capture 
into the Laplace relation as the resonance variable 9% is trapped into libration about 0° 
has a probability of about 0.9 (0). # 4 is circulating because the free eccentricity (0) of 
Ganymede's orbit exceeds the relatively small forced eccentricity in the resonance (Tjj). 
The librations are damped by the tidal dissipation in Io. Variations of this scenario have 
been proposed by Malhotra and colleagues (0, 0, |^)- The creation of the Laplace rela- 
tion by differential tidal expansion of the orbits depends on the dissipation function Qj 
of Jupiter being sufficiently small (~ 10 6 ) for non-negligible expansion of the orbits to 
occur, and the maintenance of the current configuration in equilibrium would require an 
even smaller Qj (~2x 10 4 ). 

Greenberg, motivated by the difficulty of accounting for a sufficiently small Qj nec- 
essary for the tidal origin or equilibrium described above and the observed positive hi, 
proposed a primordial origin of the Laplace relation ([|£|). He described evolution from 
much deeper in the resonances (i.e., with n\ — 2n 2 approaching zero and thereby forcing 
much larger orbital eccentricities). The path from deep resonance as dissipation decreased 
the eccentricities is characterized by stable librations of the resonance angles about centers 
much different from the 0° or 180° currently observed. The high dissipation rates at the 
high eccentricities imply rapid decay of the resonances (< 10 7 year time scale), and the 
only way to preserve them for 4.6 x 10 9 years was to have a sufficiently low Qj ~ 10 6 to 
allow at least an episodic heating and cooling of Io, albeit not as low as would be required 
by the equilibrium model. Now we would be observing the system on its way out of the 
resonance, which lowers the eccentricities and, hence, the heating rates. Eventually Io 
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cools and the dissipation function Q t of Io increases to the point at which the Jupiter 
torque "wins," pushing the system back to higher eccentricities to repeat the cycle (|3|,|i7j). 
It is not certain that Io could cool rapidly enough in this scenario, but Qj ~ 10 6 could 
delay the evolution from deep resonance sufficiently to prevent its destruction ([77]). Ma- 
jor problems with this scenario were that (i) there was no mechanism for placing the 
system in the resonances at the time of satellite formation and (ii) evolution through the 
high-eccentricity phases could not be followed with the analytic theory. Here, we describe 
a mechanism that addresses both of these problems and adds support to Greenberg's 
hypothesis. 

Constraints on the formation of the Galilean satellites of Jupiter have been established 
by theoretical analysis and the Galileo spacecraft observations. The nearly coplanar and 
nearly circular orbits of the four large satellites imply formation within a dissipative, 
equatorial disk of gas and solid particles formed, perhaps, by one of several processes 
during the formation of Jupiter flTI}). Ganymede and Callisto are about 50% ice and 50% 
rock by mass (|73|). Europa has a layer of water and ice 80 to 170 km thick ([1^). The 
large amount of relatively volatile material in Ganymede and Callisto implies that the 
disk temperature remained sufficiently cool during the accretion process that icy particles 
could persist. Callisto may not be fully differentiated (JTB, 16), a state that requires an 



accretion time scale ~ 10 6 years to keep interior temperatures below 273 K (|7^ ? [T^). 

Accretion models for the Galilean satellites have typically assumed an initial disk with 
sufficient solid material to form the satellites and with gas added to bring the mixture 
up to solar composition — the so-called minimum mass sub nebula (MMSN) (|TI|). But a 
careful evaluation of the processes in such a massive disk ( |7^ ) indicate that temperatures 
would generally be too high and accretion rates too fast to allow formation of the satellites 
with their observed properties. In particular, it would be difficult if not impossible to keep 
Callisto only partially differentiated. Therefore, it is likely that the satellites formed near 
the end of Jupiter's formation when the flux of material into Jupiter through an accretion 
disk was greatly reduced (|76|). Theoretical calculations ( p?6| ) show that Jupiter would 
limit this flux by opening a gap in the disk of material surrounding the Sun when the 
planet's mass approached its current value. Hydrodynamical simulations (|^7| ; p?I|, |^|) 
show that gas and particles continue to trickle through the gap in two streams entering 
at the inner and outer Lagrange points to form an accretion disk around Jupiter. The 
surface mass density of the accretion disk resulting from this scenario is sufficiently less 
than the MMSN to allow the slow accumulation of the satellites in a cool environment. 

By modeling processes in this disk, Canup and Ward ([71]) have found a scenario of 
satellite accretion that appears to satisfy all of the observational and theoretical con- 
straints imposed by the satellites' current properties. The important characteristics of 
the model, for our purposes, are that the formation time scales of the individual Galilean 
satellites are similar and that the satellites migrate toward Jupiter while the disk persists. 
The migration is the result of the gravitational interaction of the satellites with spiral 
density waves generated in the disk material at mean motion commensurabilities with the 
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satellites (Lindblad and corotation resonances). The satellites are not sufficiently massive 
to open gaps in the disk and, thus, they execute the so-called Type I drift Qlffi) , whose 
time scale is inversely proportional to both the disk surface mass density and the satellite 
mass. 

The key ingredient of the Canup-Ward model that allows an alternative history of the 
Galilean satellite resonant system is that Ganymede's orbit converges on those of Europa 
and Io because of its faster Type I migration induced by its greater mass. Callisto is 
left behind because of its smaller mass. Various choices of the viscosity parameter a 
and opacity of the disk lead to many different migration rates for the satellites, but all 
reasonable choices lead to converging orbits of the inner three satellites [figs. 7c, 8c, and 
9c of fll£|) l because of Ganymede's more rapid migration. Any other model of formation 
will also suffice for our purposes here, provided that the formation time scales of all the 
satellites are comparable and that the disk persists long enough for differential migration 
to lead to converging orbits 

We demonstrate (Fig. 1) that the Laplace relation can be established primordially 
during the final stages of satellite accretion by allowing Ganymede's orbit to migrate in- 
ward on the 10 5 year time scale of the Canup-Ward model fll£|), whereas Europa and Io 
initially migrate more slowly on a 2 x 10 5 year time scale. (The very rapid migration 
rates implied for disk densities corresponding to the MMSN would most likely preclude 
the capture of the satellites in any resonances during differential migration.) The satellites 
start in an initial configuration where the orbits are separated by distances greater than 
those appropriate to the 2:1 mean motion commensurability. The numerical calculations 
are done both by the Wisdom-Holman fll^ ) symplectic integrator contained in the SWIFT 
package fliffi ) and in a center-of-mass cartesian coordinate system with a Bulirsch-Stoer 
algorithm. The algorithms are modified as described in ( |I%| ) to accommodate the migra- 
tion and eccentricity damping, and both algorithms yield similar results. Europa is first 
captured into the 2:1 mean motion resonances with Ganymede and Io is subsequently 
picked up in the 2:1 resonances with Europa, simultaneously establishing the Laplace 



relation (\30j} 



The orbital eccentricities of the satellites are damped by interactions with the gas disk 
[e.g., (pIDl and eventually by tidal dissipation in the satellites [e.g., (||)]. In the case of 



disk interactions, the ratio of the inverse time scales is given by (|57| 



l/r e \e/e\ C e (H\ 

- 9tt ( — ~ 30 > (!) 



l/r a \a/a\ C a V a 

where H/a is the ratio of the scale height of the disk H to the semimajor axis of the 
orbit a. H/a « 0.1 is assumed constant in the disk, where the numerical value follows 
from the conditions in the Canup-Ward model (|T£|). C a ~ 3 is the coefficient of the 



expression for a/ a |3lD, and C e ~ 0.1 is the coefficient in the expression for e/e 
(|^|). For the integrations, we assumed that the time scales for migration in a and in 
eccentricity damping are constant and, hence, the ratio of the inverse time scales is fixed, 
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with \e/e\/\a/a\ = 30 from Eq. 1 assumed for all three satellites. The eccentricity damping 
causes the eccentricities to reach constant asymptotic values (Fig. 1), which are reduced 
if the ratio of the inverse time scales in Eq. 1 is increased. The evolution is reasonably 
independent of the rate of migration up to a certain point. An example where the satellite 
migration rates were 30 times larger led to a very similar evolution. 

Although all of the resonance variables, including the Laplace angle 9$, are librating in 
the primordial configuration of the satellites (Fig. 1), the geometry of the system differs 
from the current configuration, where all resonance variables, except #4, librate about 
either 0° or 180°, as described above. The centers of libration of 85 and 82 are near 80° 
instead of 180° and that of 8\ is near —30° instead of 0° (Fig. 2). 84 would normally 
librate about 180° for smaller eccentricities, but here it is near —10°. Only 83 has the 
same libration center near 0° that it has now. Conjunctions would no longer occur near 
the lines of apsides of the orbits. This is consistent with the locus of stable libration 
centers for 85 branching symmetrically from 180° as the system is driven deeper into the 
resonances than the current configuration ([J). The occurrence of libration centers of both 
82 and 8 5 near —80° in some of our calculations is verification of the symmetric branches. 

As the satellites move closer to Jupiter and the disk dissipates, tidal dissipation in the 
satellites eventually dominates the eccentricity damping. At the same time, the dissipation 
of tidal energy in Jupiter from tides raised by the satellites results in the torques that 
tend to push the satellites outward thereby tending to increase the eccentricities [e.g., 
([|)]. Because the current orbital eccentricities are smaller than the primordial asymptotic 
values (Fig. 1), the tendency for the eccentricities to decrease due to dissipation in the 
satellites must dominate the effect of the tides raised on Jupiter that tend to increase 
the eccentricities — at least at the end of migration. The tidal time scale for eccentricity 
damping is |e/e| ~ 6.5 x 10 4 Q I /f years for Io, which is derived from expressions given 
in ([?3p. Here Q l ~ 30 to 100 for rocky material and / is a factor to account for an 
inhomogeneous structure of Io, which can be as high as 13 for a two layer model consisting 
of a solid mantle over a molten interior (^|). The tidal time scale for eccentricity damping 
is probably more than 6 times as long as the migration time scale of the Canup-Ward 
model, and possibly 100 times as long as the disk eccentricity damping time scale if Eq. 
1 is correct. This means that tides will not take over the eccentricity damping until after 
the migration has essentially stopped. 

After migration has stopped, we assumed that the dissipation in Io is reducing its 
eccentricity on the same time scale as disk interaction did during the migration (~ 6.7x 10 3 
years) (Figs. 3 and 4). We maintained this time scale for the damping calculation to keep 
the computer time within reasonable bounds. The algorithm is such that forces are applied 
to Io that would cause the assumed eccentricity damping if Io were isolated but, because it 
must drag down the eccentricities of all the satellites locked in the resonances, the actual 
time scale is much longer. The system naturally relaxes to the current configuration 
through a series of states with offset libration centers, much like Greenberg described (fj). 
The eccentricities at the end of the calculation have almost reached the values observed 
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currently, and the libration centers have migrated to the current geometry from their 
positions at the end of migration (Fig. 2) — even to the detail that 6* 4 has ended up 
circulating just as we observe now. 

The eccentricity of Ganymede, e$, evolves in an interesting way in conjunction with 
the behavior of #4 (Figs. 3 and 4). Initially, the center of libration of 84 shifts from near 
—10° to near 180° in about 10 4 years while decreases. At 1.5 x 10 4 years, a sudden 
jump in 63, perhaps from an encounter with a secondary resonance (^), results in 84 
changing to circulation (angle increasing or decreasing without bound) where it remains 
until t = 5 X 10 4 years. Here suddenly becomes smaller, perhaps from another change 
in the effect of the secondary resonance, and #4 again librates about 180°. After ~ 10 4 
years, 84 switches from libration to permanent circulation due to the slow reduction in the 
forced eccentricity while the free eccentricity remains near 0.0016. In another calculation 
with a different migration rate, the resonance variable 84 shows similar behavior but 
without the sudden changes in and with different timing of events. In both cases, 84 
ends up circulating with e 3 ~ 0.0016. This is about the current value of Ganymede's free 
eccentricity, but there may be circumstances where less excitation of the free eccentricity 
leaves 84 librating. 

A primordial origin of the Laplace relation relieves Jupiter of providing sufficient 
tidal torque over history to expand the orbits enough for the completely tidal origin. 
However, it has been argued that even if the origin of the resonances is primordial, the 
average Qj cannot be far from the bounds 6 x 10 4 < Qj < 2 x 10 6 established by 
(\3dj , \lj\T]\,$,\$) . The current dissipation in Jupiter must still be sufficient to prevent the 
destruction of the resonances from the tidal dissipation in Io and might be sufficiently 
high to preserve the system in equilibrium (|77],[|). In contrast, could it be that we are 
witnessing the further decay of the Laplace resonance due to the dissipation in Io ([|.[|)? 
Because the energy dissipated must come from the orbits, a continuing decay of the 
Laplace configuration would result in an increasing n\. Measures of h\ — Ai, based on 
300 years of precise ephemerides of the Galilean satellites coupled with current precise 
observations, imply either an acceleration of Io's mean motion, indicating further decay of 



the resonant configuration (|^,P3), or insufficient deceleration (|5^) from Jupiter torques on 



Io to accommodate the high heat flux from Io in an equilibrium eccentricity configuration. 
The heat flux is now estimated to be about 3W/m 2 §40j). Although the measurements 
of hi remain controversial, we should keep in mind that a heat flux from Io of 3W/m 2 
would require a current Qj ~ 2 x 10 4 to maintain the system in equilibrium, which is a 
factor of 3 lower than the minimum historical average Qj = 6 x 10 4 from the proximity 
of Io to Jupiter (|ff6|) . But there is no reason to believe that the current Qj should be 
constrained much by a bound on its historical average (Q). 

A primordial origin of the Laplace relation among the Galilean satellites, created by 
differential migration of the satellites in the gaseous disk, with the help of tidal dissipation, 
reproduces every detail of the current configuration. Changing assumptions will change 
details of this scenario such as the asymptotic values of the eccentricities during the 
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migration; a different damping rate may change the circulation of 64. But as long as the 
satellites exist simultaneously in the disk and the more rapid migration of Ganymede forces 
the inner-three orbits to converge and capture into the Laplace resonance results 
the qualitative features of this alternative early history will persist. 

References and Notes 

1. C. F. Yoder, Nature 279, 767 (1979). 

2. C. F. Yoder, S. J. Peale, Icarus 47, 1 (1981). 

3. R. Greenberg, in Satellites of Jupiter, D. Morrison, Ed. (Univ. Arizona Press, Tucson, 
AZ, 1982), pp.65-92. 

4. R. Greenberg, Icarus 70, 334 (1987). 

5. Free and forced eccentricities are often referred to in discussions of eccentricity type 
orbital resonances at commensurabilities of the mean motions. In a lowest order ec- 
centricity type resonance between two satellites, where for example the resonance 
variable Ai — 2\ 2 + z&i is librating about 0° {ni — 2n 2 + w\ = on average), conjunc- 
tions must occur when satellite 1 is near its periapse for stability. This imposes the 
constraint (wi) = — (n x — 2n 2 ), where the angled brackets denote averaged values. 
Because {w\) depends on the magnitude of the eccentricity e\, the resonant interac- 
tion will force e\ to the appropriate value to accommodate the resonance. If there is 
a nonzero libration amplitude, the eccentricity will fluctuate about the forced mean 
value. The magnitude of the fluctuation in e± is called the free eccentricity, which can 
be damped to zero. 

6. S. J. Peale, Annu. Rev. Astron. Astrophys. 37, 533 (1999). 

7. R. Malhotra, Icarus 94, 399 (1991). 

8. A. P. Showman, R. Malhotra, Icarus 127, 93 (1997). 

9. A. P. Showman, D. J. Stevenson, R. Malhotra, Icarus 129, 367 (1997). 

10. G. W. Ojakangas, D. J. Stevenson, Icarus 66, 341 (1986). 

11. S. J. Peale, R. Greenberg, Lunar Planet. Set. XI, 871 (1980). 

12. J. B. Pollack, J. I. Lunine, W. C. Tittemore, in Uranus, J. T. Bergstralh, E. D. Miner, 
M. S. Matthews, Eds. (Univ. Arizona Press, Tucson, AZ, 1991), pp. 469-512. 

13. G. Schubert, T. Spohn, R. T. Reynolds, in Satellites, J. A. Burns, M. S. Matthews, 
Eds. (Univ. Arizona Press, Tucson, AZ, 1986), pp. 224-292. 



7 



14. J. D. Anderson et al, Science 281, 2019 (1998). 

15. J. D. Anderson et al, Science 280, 1573 (1998). 

16. J. D. Anderson et al, Icarus 153, 157 (2001). 

17. D. J. Stevenson, A. W. Harris, J. I. Lunine, in Satellites, J. A. Burns, M. S. Matthews, 
Eds. (Univ. Arizona Press, Tucson, AZ, 1986), pp. 39-88. 

18. D. J. Stevenson, Science 294 71 (2001). 

19. R. M. Canup, W. R. Ward, Astron. J., in press. 

20. W. R. Ward, J. M. Hahn, in Protostars and Planets IV, V. Mannings, A. P. Boss, S. 
S. Russell, Eds. (Univ. Arizona Press, Tucson, AZ, 2000), pp. 1135-1158. 

21. S. H. Lubow, M. Seibert, P. Artymowicz, Astrophys. J. 526, 1001 (1999). 

22. G. Bryden, X. Chen, D. N. C. Lin, R. P. Nelson, J. C. B. Papaloizou, Astrophys. J. 
514, 344 (1999). 

23. W. Kley, G. D'Angelo, T. Henning, Astrophys. J. 547, 457 (2001). 

24. G. D'Angelo, T. Henning, W. Kley, Astron. Astrophys. 385, 647 (2002). 

25. W. R. Ward, Icarus 126, 261 (1997). 

26. Lee and Peale (^) have shown that capture into the resonances at the 2:1 commen- 
surability during differential migration of the two planets in coplanar orbits about the 
star GJ 876 occurs easily for a reasonably wide range of initial conditions. Capture 
and subsequent configuration are independent of the migration rate unless that rate 
is unreasonably large, and the capture and configuration are independent of the func- 
tional dependence of the migration on orbital semimajor axes. The only thing that 
matters is that the orbits converge. 

27. J. Wisdom, M. Holman, Astron. J. 102, 1528 (1991). 

28. H. F. Levison, M. J. Duncan, Icarus 108, 18 (1994). Available at: 
www . b oulder . swri . edu/ ~ hal/swif t . ht ml . 

29. M. H. Lee, S. J. Peale, Astrophys. J. 567, 596 (2002). 

30. We have arbitrarily assigned the full masses to the satellites and assumed that Io 
and Europa migrate on the same time scale to demonstrate the hypothesis, whereas 
migration will ensue as the satellites are accreting and Io will generally migrate faster 
than Europa once its mass exceeds that of the latter. The necessary condition here 



8 



for the success of the establishment of the Laplace relation is that Io migrate more 
slowly than Europa after Europa is captured into resonance with Ganymede, which 
capture accelerates Europa's migration. The migration rate history thus depends on 
the surface mass density of the disk and the accretion history determining the satellite 
masses as a function of time. There is enough flexibility in the disk and accretion model 
to effect the necessary relative rates. 

31. P. Artymowicz, Astrophys. J. 419, 166 (1993). 

32. H. Tanaka, T. Takeuchi, W. R. Ward, Astrophys. J. 565, 1257 (2002). 

33. S. J. Peale, P. Cassen, R. T. Reynolds, Icarus 43, 65 (1980). 

34. S. J. Peale, P. Cassen, R. T. Reynolds, Science 203, 892 (1979). 

35. W. C. Tittemore, J. Wisdom, Icarus 78, 63 (1989). 

36. P. Goldreich, S. Soter, Icarus 5, 375 (1966). 

37. S. J. Goldstein, K. C. Jacobs, Astron. J. 110, 3054 (1995). 

38. K. Aksnes, F. A. Franklin, Astron. J. 122, 2734 (2001). 

39. J. H. Lieske, Astron. Astrophys. 176, 146 (1987). 

40. G. J. Veeder, D. L. Matson, T. V. Johnson, A. G. Davies, D. L. Blaney, AGU Spring 
Meeting 2002, 28 to 31 May 2002, Washington, DC (American Geophysical Union, 
Washington, DC, 2002). Available at: www.agu.org/dbasetop.html. 

41. The establishment of the Laplace relation is not a certain outcome for differential 
migration of the satellites, although some type of resonance configuration does seem 
certain. Minor changes in the initial conditions, migration rates, or eccentricity damp- 
ing rates sometimes led to other configurations. For example, in one case the system 
ended up with Ganymede-Europa in 2:1 resonances and Europa-Io in 7:3 resonances. 
Accretion of the satellites may not have been complete at the time of capture into the 
resonances, but two identical calculations — one with the full satellite masses and the 
other with half these masses — led to very similar Laplace configurations. Even if the 
mass ratios vary during migration, the configuration at the end of migration should 
be the same. The conditions for the various end configurations and an estimate of 
the probability of the system selecting the Laplace configuration will be subjects of a 
future paper. 

42. We thank R. Canup and W. Ward for sending us an early preprint of their paper on 
Galilean satellite origin. This research was supported in part by NASA PGG grants 
NAG5 3646 and NAG 5 11666 and SSO grants NAG5 7177 and NAG5 12087. 



9 




5xl0 4 10 5 

t (years) 



Figure 1: Nebula-induced evolution of Galilean satellites into the Laplace resonance. The 
semimajor axes are normalized by Ganymede's current distance from Jupiter 0,3^. The 
imposed inward migration and eccentricity damping time scales are 1 03/03! ~ 10 5 years, 
= I02/02I ~ 2 x 10 5 years, and | / | = |aj/(30dj)| % = 1, 2, and 3. Circular 
coplanar orbits are assumed initially. 
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Figure 2: Libration of 2:1 mean motion resonance variables and the Laplace angle for the 
asymptotic state of the resonant Galilean satellite configuration shown in Fig. 1. The 0j 
are defined in the text. The coefficients ej corresponding to the Oi are indicated in the 
figure. Librations are shown for 250 years with all migration and damping removed. 
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Figure 3: Relaxation of Galilean satellite system to the current configuration with mi- 
gration halted and with eccentricity damping in Io alone to simulate the effects of tidal 
dissipation. 
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Figure 4: Behavior of resonance variables for the evolution described in Fig. 3. 
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